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Leptonic D and Ds Decays near cc Threshold 

S. R. Blusk 

Syracuse University, Syracuse, NY 13244, USA 



We present recent results from the CLEO Collaboration on leptonic decay rates of D and Dg near cc production 
threshold. From these decay rates, we extract the decay constants, /q+ = (222.6 ± 16.7_^'^) MeV, f + = 

(274 ± 10 ± 5) MeV, and the ratio fj^+/f]j+ = 1.23 ± 0.11 ± 0.03. 
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1. Introduction 

Within the Standard Model, leptonic D (or B) me- 
son decays proceed via annihilation of the initial state 
quarks. The matrix element is described by the prod- 
uct of a hadronic current, a leptonic current, along 
with a W propagator. The form of the latter two 
are well-known within the Standard Model, however, 
the hadronic matrix element, which represents the an- 
nihilation of the initial state heavy quark and light 
antiquark, depends on the details of the initial-state 
quark wave- functions, and is not calculable using stan- 
dard techniques of perturbative QCD. This hadronic 
matrix clement can be computed using either lattice 
QCD 1, 2, 3], or other techniques 0, H i, 0, i • The 
partial width for the leptonic decay is given by: 
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is the D+ 



where Gp is the Fermi constant, 

(Df) mass, mi is the final state lepton's mass, and 
Vcd(s) are the relevant CKM matrix elements. The 
quantity fu is the decay constant and represents the 
hadronic matrix element discussed above. A critical 
input to B mixing and CP violation measurements in 
the B sector is the B decay constant, Jb- Due to 
the difficulty in measuring /b , we take the value from 
theory, usually lattice QCD. To have confidence in 
the theoretical number, a stringent theoretical test is 
provided by a precision measurement of the D decay 
constant, fo- Such a measurement provides a critical 
test of any theory or model that makes predictions for 
decay constants. 

The CLEO experiment, operating near cc threshold, 
is well positioned to measure these decay rates, and 
hence /£,+ and /£>+• Charge conjugate finals states 
are implied throughout unless otherwise noted. 



the production threshold implies that the ip{377Q) de- 
cays to DD with no additional particles. We exploit 
this clean final state, along with the hermiticity of the 
detector to reconstruct the neutrino from the miss- 
ing momentum in the event. Specifically, we fully re- 
construct a D^ meson (the tag) in six hadronic fi- 
nal states, comprising A^tag = 158, 354 ± 496 tags. 
To search for D'^ — > /i^J^, we require a single ex- 
tra charged particle with an energy deposition in the 
crystal calorimeter (CC), E^^ <300 MeV, and veto 
events with any additional photon candidates with 
energy larger than 250 MeV. From this subsample 
of events, we compute the square of the missing- 
mass {AdM'^) recoiling against the D~ jX^ system. For 
D^ — > ^i'^h'fi, a peak at zero is obtained with a res- 
olution of a{MNP) - 0.025 GeV^. The MM"^ dis- 
tribution is shown in Fig. [T] for data. The clear 
excess near zero is the D^ -^ ji^v signal. Some 
D"*" -^ Ks^lt^^ events pass the selection requirements 
and appear as a prominent, but well-separated peak 
near MM'^ ~ 0.25 GeV^. 
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2. Measurement of / 



D 



To measure /^ Q, we use 281 pb ^ of data col- 
lected at the -0(3770) resonance. The proximity to 



Figure 1: Missing-mass squared distribution for D — > 
jj.'^v candidates. The peak near zero corresponds to signal 
events, and is expanded in the inset. The larger peak at 
MM'^ ~ 0.25 GeV^ corresponds to D+ -^ Ks,Ln'^ events 
which pass the selection requirements. 
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The branching fraction is computed using: 
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where iVcand = 50 is the number of signal candi- 
dates in the region \MM^\ < 0.050 GeV^ iVback = 
2.81 ± 0.30 ± 0.27 is the expected number of back- 
ground events, iVtag = 158, 354 ± 496 is the number 
of fuUy-reconstructed D~ tags, e^ = 69.4% is the ef- 
ficiency for reconstructing and identifying the muon, 
and ecc = 96.1% is the fraction of events that do 
not have any additional photon candidates with en- 
ergy larger than 250 MeV. An additional correction of 
(1.5 ± 0.4 ± 0.5)% is applied to account for the higher 
efficiency for reconstructing a D~ tag in D+ —>■ ^^v 
events than in generic hadronic events. 
The resulting branching fraction is 



B{D+ -^ ^i+v) = (4.40 ± 0.66J:o:?2) x 10"^- (3) 
Using Eq. [1] we determine the decay constant to be: 
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(222.6 ± 16.7tl:l) MeV. 



(4) 



3. Measurement of / 



d: 



The measurements of /^+ at CLEO require higher 
energy running in order to produce the DgDs pair. 
A scan of the energy region from 3970 to 4260 MeV 
was performed, and it was determined that the opti- 
mal energy for Ds physics was 4170 MeV [11|, where 
D,D* is dominant, e.g., a{D,Dl) = (916 ± 50) pb 
and a{DsD*) = (35 ± 19) pb. A slight complication 
with using DgD* is the additional (--^150 MeV) pho- 
ton(s) from the D* decay. Two independent analyses 
have been carried out. The first analysis is similar 
to the Z?+ -^ 11^ ly measurement described previously, 
where, in addition to measuring B{D~^ —>■ fJ.'^i'), we 
also measure B{Df — > t^v), where, t+ -^ tt^vv. 
In the second analysis, we measure B{Df -^ t^i^), 



3.1. Measurement of i3(D+ -^ {fi^ 
using Missing Mass 



.T^ ly 



We use 314 pb^^ of data collected at Ecm = 
4170 MeV for this analysis. We search for final states 
consistent with either Df — > /x+i' or Df -^ t'^v. The 
branching fraction is obtained from: 
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where N'^ is the number of reconstructed DgD* 
events and e is the efficiency for reconstruction and 
identification of the /i''" for D+ — > /i"*"!^, or the 7r+ for 
Df -^ T^v, r+ -^ TT^vrTu We therefore absorb the 
full reconstruction of the D* into the denominator, 
and do not rely on Monte Carlo simulation for the 
efficiency of the ~100 MeV photon. 

To determine N^ , we first fully reconstruct a 
hadronic Dj tag in eight tag modes, from which we 
obtain 31,302 ± 472 D+ tags. To identify DsD*+s 
events, we combine a D~ tag with any additional pho- 
ton candidate in the event and form the missing-mass 
squared {MM*"^) recoiling against the ^D'^ system, 
MM*2 = [E,^ - Ed^ - E.,f - (pcm - PD, - Pj?. 
This quantity peaks at Af|, , regardless of whether 
the photon came from the Dj (the tag) or from the 
Df. The distribution of MM*"^ is shown in Fig [2] for 
all eight tag modes combined. A fit to this distribu- 
tion yields 18645 ± 426 ± 1081 0^0^ events within 
±2.5 standard deviations of M{Dg). 
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Figure 2: Square of the missing mass recoiling against a 
7D* candidates. 



To search for ZJj 



^i^v 



and D+ 



'K^vnu, we require a single additional charged parti- 
cle and no additional photon candidates with energy 
in excess of 250 MeV. The signatures for D+ -^ [i^v 
and Z)+ — > T'^v, t+ — > -n^vnu are similar in that 
they both have a Dg tag and a single high momentum 
charged particle. In addition to the difference in the 
energy depositions of muons and pions, the two-body 
versus three-body decay implies significantly different 
missing mass (MM"^) distributions. To suppress back- 
grounds with neutrals, we veto events which have an 
energy deposition (excluding the tag) in the CC ex- 
ceeding 250 MeV. The two-body leptonic decay form a 
MM^ distribution that peaks near zero with a resolu- 
tion of ~0.025 GeV^. The three-body leptonic decay 
covers a broad MM"^ region, which peaks near 0.1 
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GeV^, and falls smoothly to zero at MRP = -0.05 
and extends to MAP -- 0.8 GeV^ We thus define 
signal samples as follows. (i)-/i: For D+ — > /x+i 
require an energy deposition, E^^ < 300 MeV 



, we 
and 



\MM^\ < 0.05 GeV2. For D+ . , ., , . „ .., 
we define two subsamples - (i)-T: E^^ < 300 MeV 
and 0.05 < MIVP < 0.20 GeV^, and (ii)-r ^^'^ > 
300 MeV and -0.05 < MAP < 0.20 GeV^. The 
upper cutoff in A'lA'P is to avoid background from 
Dt -^ K^TT^ . We also consider a third sample, (iii)- 



e, 



for D+ 



e^v by requiring the track's energy 



deposition to be consistent with its momentum and 
|MM2| < 0.050 GeV^. 

The AIA'P distributions are shown in Fig. [31 where 
cases (i)-/.t and (i)-T are combined. In the Df —i- 
H+v signal region of jAfM^I < 0.05 GeV^ we find 
92 events with an expected background of 3.5 ± 1.4 
events. This sample is mostly of Z)+ -^ ^-^v, with 
some cross- feed from _D+ -^ t^v, r+ — > ti^vD. We 
evaluate the branching fraction using the number of 
/i+z/ events, N^^^92, in the signal region: 



N,,,, = 
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where e — 80.1% is the efficiency of reconstructing 
the charged particle in a D^ — > ^.^v event, and in- 
cludes the veto on events with additional photons with 
E > 250 MeV. The quantity, e' = 91.4%, is the prod- 
uct of the muon identification efficiency (99%) and 
the MM^ < 0.05 GeV^ requirement (92.3%). The 
cross-feed efficiency, e"=7.9%, which is the product 
of the efficiency of the pion depositing less than 300 
MeV in the CC (60%) and the MAP < 0.05 GeV^ 
requirement (13.2%). One can re-express B{D^ -^ 
T'^v, T+ — * Ti^vD) as: 



20 


1770107-011 

(i) 


15 




; 


10 


■ F 


i-r 






■ 


5 

CM 
> 


ir 1 


J u 


knjln ^, 

, 1 , 1 , , , , 


i" 


. (ii) 


w 




c 


■ 




> 




LU 
5 


i n n r 






1 
n 


-(iii) nn M n|-[ 



0.25 

MM^ (GeV^) 



0.50 



Figure 3: Square of the missing mass recoiling against 
7D*/i+(or 7r+) candidates for cases (i) Eqq < 300 MeV 
(_D^ — > /i^i/ and D^ -^ t'^v, t^ — > ■k'^uv candidates com- 
bined), (ii) Ech > 300 MeV (D+ -^ r+v, r+ ^ -k+vv), 
and (iii) charged particle consistent with an electron. 



B{Dt 



Ti+vv) = R ■ i3(T^ 



xB{Dt -^ ii+v) = 1.059 • B{Dt -^ ^i'^v) (7) 
where we use the Standard Model ratio for R: 



R 



not 



1- 



TiD+ -^ ^i+v) 



mT-+ 



1- 



We thus find: 



9.72. 



(8) 



B{Dt -^ i-i+iy) = (0.594 ± 0.066 ± 0.031)%, (9) 

where the 5.2% systematic error is dominated by the 
5% uncertainty on N^^^. 



We also compute B{Df -^ r+t/, t+ -^ tt+i^j/) us- 
ing cases (i)-T and (ii)-T. For these two cases, we 
find yields of 31 and 25 events, and expected back- 
grounds of 3.5l|^'^ and 5.1±1.6 events, respectively. 
The fraction of Df -^ t^v, r+ — > -k'^vv events in the 
respective MAP regions are 32% and 45%. We thus 
find: 



B{Dt 



vv) 



(8.0 ±1.3 ±0.4)%. 
(10) 
With the measured branching fractions, B{D^ -^ 
/x+z/) and B{Df —>■ t~^v, r"*" —f tt'^vv), we measure 
the ratio of partial widths, R=13.4±2.6±0.2 (defined 
in Eq.[5]), which is consistent with the Standard Model 
value of 9.72. 

We may improve on the precision of B{Df -^ fJ'^i') 
by combining the Df — > fi'^v and Df — *■ t'^v, t+ -^ 
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TT+z/i/ candidates. We can still use Eq. [71 except e' 
and e" increase from 91.4% and 7.9% to 96.2% and 
45.2%, respectively. We thus find an effective branch- 
ing fraction: 



= 1.23 ±0.11 ±0.04. 



ID+ 



(13) 



B^'^iDt -^ ^x+v) = (0.638 ± 0.059 ± 0.033)%. (11) 4. Measurement of D+ -^ t+u, t+ -^ e+vv 



Again, the dominant systematic uncertainty (5%) is 
on the number of D* tags. 

The MAP distribution for all selected Df -^ ^l^v 
and D^ —> t^u, t"*" — > ■n'^vv candidates is shown 
in Fig. m Overlayed is a curve that represents the 
expected shape, normalized to the event yield in the 
data in the MM'^ region below 0.2 GeV^. We find 
good agreement between the shape in data and ex- 
pectations. 
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In the second measurement of B{Df -^ t^v), we 
use 298 pb"^ of data collected at E^^ = 4170 MeV. 
We utilize the decay t"*" — > e'^vv^ where we bene- 
fit from the large value of ,B(t+ -^ e'^vv) ~ 18%, 
and the excellent electron identification capabilities of 
the CLEO-c detector. We fully reconstruct the three 
hadronic decay channels: Dj -^ (j)Tr~ ,K*^K~ and 
KgK~. Charged hadrons are identified using stan- 
dard selection criteria [12l |. and the intermediate res- 
onances, (j) -^ K+K-, K*° -^ K-TT+, and K^ -^ 
TT~^Tr~ , are required to have an invariant mass within 
±10 MeV, ±75 MeV and ±12 MeV of their known 
values [131 . Signal candidates are required to a recon- 
structed invariant mass, M{Ds) within ±20 MeV of 
the known Dg mass (tod^, ). We also define sideband 
regions, 35 < |M (£)«) - mo, \ < 55 MeV, to study the 
combinatorial background. The invariant mass distri- 
butions of the three D^ tag channels are shown in 
Fig. El 
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Figure 5; Invariant mass distributions of DJ candidates 
from data. The points are data, the solid line is a fit, and 
the dashed line is the background. 



Figure 4: Square of the missing mass recoiling against 
7Ds/^'''(or TT^) candidates. The curve is the expected 
shape from simulation, normalized to the number of events 
with MM'^ < 0.2 GeV. 



We also search for the decay -D+ 



The he- 



licity suppression in this decay is much larger, and 
the expected rate is ~50,000 times smaller than in 

v candidates and 



-^ e^ 
e'^iy 



< 1.3 X 10- 



at 



D+ -^ jjL+v. We find no D+ 
set the upper hmit, B{Df — 
the 90% confidence level. 

Using the more precise value for B{D^ -^ l-i-^'^) 
from Eq. [Til we compute the decay constant, fj^+y. 



fr 



274 ± 13 ± 7 MeV 



(12) 



Combining this with our previous result for fjj+ = 
(222.6 ± 16.713-4) MeV we determine the ratio: 



To ensure we have D^D*, we compute the mass re- 
coiling against the reconstructed Dg , and require it to 
be within ±55 MeV of the D* mass ^. We then 
select the subset of events with a single additional 
charged track with p > 200 MeV that has opposite 
charge to the Ds tag and is consistent with being a 
positron. The discriminating variable we use to iden- 
tify Df -^ T^i/, T+ -^ e^vV is -Eextraj the total en- 
ergy remaining in the calorimeter after all showers as- 
sociated with the tag and the positron are removed. 
In signal events, the only additional particles beyond 
the Ds tag and the positron are the two neutrinos 
and either a photon from D* -^ 7-Ds, or a tt*^ from 
_D* -^ TT^Ds ■ Kinematically, these photons populate 
the energy regions from 114-170 MeV (for 7-D+) and 
39-117 MeV {hounT^D+). 

The distribution of Soxtra in data is shown in 
Fig. [HI The large excess at low values of -Ecxtra 
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is the £>+ -^ T^v, t+ -^ e^vv signal. The 
broad background which peaks near 1 GeV is pre- 
dominantly semi-leptonic decays, such as Dj" -^ 
^e^v, r]e+v, r(e^v, K^e+v and K*^e+v. The 
Cabibbo-suppressed decay, K^e'^v, produces a small 
peaking component in the signal region. The shape 
of this background is taken from Monte Carlo sim- 
ulation, and is normalized to our measured rate for 
D+ ^ Kle+v of B{Dt -^ Ks^'^iy) = (0.27 ± 0.10)%. 
We choose the signal region as -Eoxtra < 400 MeV, 
which is chosen based on optimizing the signal signif- 
icance. The expected non-peaking background in the 
signal region is estimated by scaling the number of 
data events with £^oxta > 600 MeV by the MC ratio of 
events in the sideband (-E^^a > 600 MeV) to signal 
region {E^^^e. < 400 MeV). The yields of Dj tags and 



Dt 



T+Z/, 



e"*";^;^ signal events are shown in 



Table m The scale factor, s shown in Table |I] is a cor- 
rection to account for slight differences in the expected 
number of background events in the signal and side- 
band regions. Using the efficiency to reconstruct the 



final state, Df 
and the B{t+ - 



T~^V, 



of (71.4 ± 0.4)% 



>) = (17.84 ±0.05)%, we find: 



B(D+ -^ Tv) = (6.24 ± 0.71 ± 0.36)%. (14) 

The 5.8% systematic uncertainty is dominated by the 
4.3% contribution fron 
ing in the calorimeter. 



4.3% contribution from the simulation of K^ shower- 
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Figure 6: Total extra energy left over in the calorimeter 
after removing energy associated with the Ds tag and the 
positron. Data are shown as the points with error bars, 
and the MC background predictions are shown as solid, 
dashed and hatched histograms, and the expected signal 
contribution is indicated by the shaded histogram. 

Using Eq. [fl we find /^+ = (275 ± 16 ± 8) MeV. 



When this result is combined with the result in Eq. [T2l 
we obtain: 



>D2 



274 ± 10 ± 5 MeV 



(15) 



5. Summary 

We have presented measurements of the branch- 
ing fractions B{D^ — > M^'^i I^t '^ M^^ ^-nd D^ -^ 
T^v, r"*" -^ TT^vv with the CLEO-c detector. The 
results are the most precise measurements of these 
leptonic decay rates to date. Using EqlU we extract 
the decay constants: 



fD+/fD+ 



(222.6 ± 16. 7l3^) MeV. (16) 
(274 ± 10 ± 5) MeV (17) 

1.23 ±0.11 ±0.03 (18) 



Our measurement of /^+ is consistent with and sig- 
nificantly more precise than the recent measurement 
by BaBar [lij. The only other measurement of /£>+ 
was reported by BES based on 1 signal candidate. Re- 
cent lattice QCD predictions [iaj[3 of both /£,+ and 
fjj+ are typically ~10% lower than our measurements, 

whereas the ratio of f]^+/fD+ is in good agreement 
with our measurement. 

We gratefully acknowledge the effort of the CESR 
staff in providing us with excellent luminosity and 
running conditions. We also thank the National 
Science Foundation for support of this research. 
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Table I Summary of D^ tagged events (yield, background from sidebands, sidebands scale factor (s), and sideband- 



subtracted yield), and Z)+ -^ r^v, r+ — + e^vv events (yield, background from D^ 
semileptonic decays, and sideband-subtracted yield). 



sidebands, background from D^ 



Mode 


Ds Tags 


D+ -^ T^v, r+ -^ e'^uv 


Yield 


Back 


s 


Signal 


Yield 


Background 


Signal 


DJ 


Dt 


D~ ^ (t>Ti- 
D- -^ K-K*° 
D- ^K-Kl 


5232 
8937 
3468 


388 
3618 
695 


1.001 
1.008 
1.030 


4843.6 ± 75.0 
5289.2 ±112.2 
2751.8 ±64.7 


49 
55 
24 



3 
2 


8.7 ±0.6 
8.5 ±0.7 

3.8 ±0.4 


40.3 ±7.0 
43.5 ±7.6 
18.1 ±5.1 


Total 


17637 


4701 


- 


12884.6 ± 149.7 


128 


5 


21.0 ±1.0 


101.9 ±11.5 



